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In this talk we discuss how the down type left-right squark mixing in Super- 
symmetry can induce a new source of CP violation in the time dependent 
asymmtries in B — > <f>K process. We use QCD improved factorization pro- 
cess to calculate the hadronic matrix element for the process and find the 
allowed parameter space for p and <p, the magnitude and phase of the down 
type LR(RL) squark mixing parameter S b ^ R , RL y In the same allowed regin we 



calculate the expected CP asymmtries in the B — > <j>K* process. 
1 Introduction 

Time dependent asymmetries measured in the decay B — > <pKg both by BaBar 
and Belle collaborations [21 ED E] show significant deviation from the stan- 
dard model and this has generated much theoretical speculation regarding 
physics beyond the standard model |5I21| . In the standard model, the process 
B — > 4>Ks is purely penguin dominated and the leading contribution has no 
weak phase. The coefficient of sin( Amst) in the asymmetry therefore should 
measure sin 2(3, the same quantity that is involved in B — > ipKs in the stan- 
dard model. The most recent measured average values of asymmetries are 



The value for S^k s agrees with theoretical expectation from the CKM matrix 
of S^Ks = sin_2/3 = 0.715to °4 5 |2IJ This leads to the conclusion that CP 
phase in B — B mixing is consistent with the standard model. The deviation 
in the 4>Kg is intriguing because a penguin process being a loop induced pro- 
cess is particularly sensitive to new physics which can manifest itself in a loop 
diagram through exchange of heavy particles. In this talk |2U we consider 
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= 0.734 ±0.055 
S (t>Ks = -0.15 ±0.33 
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effects arising from non universal squark mixing in the second and third gen- 
eration of the down type squarks in supersymmetric theory as the origin of 
additional contributions to the amplitude within the mass insertion approxi- 
mation scheme. In particular, the exchanges of gluinos (g) and squark (q) with 
left-right mixing can enhance the Wilson coefficient of the gluonic dipole pen- 
guin operator Os g by a factor of m§ / m& compared with the standard model 
prediction and we take into account its effect on the process B — > <j>Ks- In 
our analysis we take the B — B mixing phase the same as in the standard 
model as required by tpKs data, and permitted in SUSY by requiring that 
the first and third generation squark mixing to be small. We study B — ► (pK 
in QCD improved factorization scheme (BBNS approach ) |25|. This method 
incorporates elements of naive factorization approach (as its leading term ) 
and perturbative QCD corrections (as sub-leading contributions) and allows 
one to compute systematic radiative corrections to the naive factorization for 
the hadronic B decays. 

In supersymmetry, assuming masses of squarks (q) and gluinos (g), the new 
source of CP violation can be parameterized by the complex quantity 5^ R ^ RL ^ 
written in the form pe 1 ^ We identify the region in p — ip plane allowed by the 
experimental data on B — > <j)K time dependent asymmetries S^Ks an d C^Ks 
and the branching ratio. This allowed region is dependent of the QCD scale //, 
therefore we illustrate the region for two values of p = m,b and mf,/2. The same 
contribution should also be present in other penguin mediated process. We 
study the effect of LR(RL) mass insertion to the B —> cj)K* decay mode which 
is also a pure penguin process using QCD improved factorization method. We 
then estimate the branching ratio B(B — > <pK*) and the CP asymmetry Acp 
in the parameter space of &y!b.(rl) a U° we d by B — > <f>K data. In this vector 
vector final state, one can also construct more CP violating observables |26| . 
We compute these observables in the same range of parameter space as that 
allowed by B — > cj)K . 



2 CP Asymmetry of B — > <pK 

The time dependent CP asymmetry of B — > <pKs is described by : 

= giMzMh^) (2) 

r(B°(t) -► 4>K S ) + r(B°(t) -► cf>K s ) 
= -C<j,Ks cos(Am B t) + S^k sm(Am B t) (3) 

where Sqk and C^k s are given by 

21m X^Ks n 1- I ^4>k s I 2 , A s 

s ' k -i+\x,k s r * Ks ~ i+ 1 hK S i 2 (4) 

and X,pKs can be expressed in terms of decay amplitudes: 
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X<t>Ks ~ " e M{W^Ws) { ) 

The branching ratio and the direct CP asymmetries of both the charged 
and neutral modes of B — ► <j>K have been measured [Tl 121 131 111 1221 177) : 

B(B° 
B(B+ - 



A C p(B^ 



3 The exclusive B — > <f>K decay 

In the standard model, the effective Hamiltonian for charmless B — ► <j)K{<f>K*) 
decay is given by 

C \ 10 

W e// = —LVa,VZ{ + tf 2 (/i)0 2 (/i) + £ WW + C 7l 7l 



<j>K s ) 


= (8.0 ± 1.3) x 1CT 6 


(6) 


4>K+) 


= (9.4 ±0.9) x 10~ 6 , 


(7) 


S<pK s 


= +0.45 ± 0.43 ± 0.07 (BaBar); 


(8) 




= -0.96 ±0.50+°;°? (Belle); 


(9) 


C^K S 


= -0.19 ±0.30 


(10) 


<f>K+) 


= (3.9 ±8.8 ±1.1)% 


(11) 




i=3 



(12) 



where the Wilson coefficients Cj(/z) are obtained from the weak scale down 
to scale fi by running the rcnormalization group equations. The definitions of 
the operators and different Wilson coefficients can be found in Ref. |25| . 



4 B -> cf)K in the QCDF Approach 



In the QCD improved factorization scheme, the B — > cpK decay amplitude 
due to a particular operator can be represented in following form : 

<<t>K\0\B >=< <t>K\0\B> fact + r " a " + °( A QCD/m b )\ (13) 

where < <fiK \ O | B >f ac t denotes the naive factorization result. The second 
and third term in the bracket represent higher order a s and AQCD/m-b correc- 
tion to the hadronic transition amplitude. Following the scheme and notations 
presented in Ref. |28ll2"9"] . we write down the B — * (f>K amplitude in the heavy 
quark limit. 
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M(B+ -> 4>K+) = 

M(B° <^°) = 9± m l UF B^K [m T )VpbV;s [a p +< + a P 

K + a 9+ a lo) | n P 

2 a 10a 

where p is summed over w and c. The coefficients af are given by 



a% = a 3 ] = C 3 + ^- 
G 3 



1 + + 
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a£ = a§ = G 5 + ^| 
a« = a c 7 = C 7 + ^ 
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Q<p 



(14) 



(15) 



with C F = (7V C 2 - l)/2N c and N c = 3. The quantities Vj,, H^, P$ and are 
hadronic parameters that contain all nonperturbative dynamics. 



JL 

m b 



^ = -121n^-18 + 4, 



/■! 1 — 2x 

fi= dxg{x)^{x)\ g(x)=3- In 

Jo ix 



x — 2>m, 



Ha, = 



4^ 2 Sb!k f 1 dz $B(z) f 1 ^k(x) f 1 



N c F»-> K (0)m% Jo 



dx- 
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Pl = G 3 [G^ 8a ) + Gt(s b )} + C 2 G4s p ) + (G 4 + C 6 ) &4>(*f) + Cf g f G 4 

f=u 

3 6 3 



G (s) = GM (2/3) 
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dx 



(l-x) 



du u (1 — u) In [s — u(l — u)(l — x)] 



(16) 



New Source of CP violation in B physics ? 



5 



where, Sj = m 2 /m 2 . Here, represent contributions from the vertex cor- 
rection and correspond to hard gluon-exchange interactions with spec- 
tator quarks. and represent QCD penguin contributions. We neglect 
order a em EW penguin corrections to aj. fs,fK are the B and K meson 
decay constants and F^^ K denotes the form factor for B — > K transitions. 
<Pb(z),<Pk(x), and &<f>(y) are the B,K, and <f> meson wave functions respec- 
tively. In this analysis we take following forms for them [2E] 

$ B {x) = N B x 2 (1 - xf exp 

*kA x ) = 6 x (1 - a;) (17) 

where, Nb is a normalization factor satisfying dx <L> B (x) — 1, and lo b = 0.4. 

For the sake of completeness, we give the branching ratio for B — ► <f>K 
decay channel in the rest frame of the B meson. 

BK(B -> cj)K) = | M(B -> cj>K) | 2 (18) 

where, tb represents the B meson lifetime and the kinematical factor | P cm | 
is written as 



I Pcm \ = 7, V [ m B - ( m K + m 4>) 2 ] \ m B ~ ( m K ~ "ty) 2 ] (19) 



2u% 



5 SUSY gluino contributions to B — > 4>K 

In order to study the new physics contribution to the CP violating phase of 
amplitude M. (B — > <f>K) , we compute the effect of flavor changing contribution 
to B — > (f>K arising from q — q — g interactions in supersymmetric theory under 
the mass insertion approximation scheme |301 131j . In this approximation, the 
flavor changing contribution is parameterized in terms of 8 l ^ B = A z l B /fh , 
where, A represents the off-diagonal entries of the squark mass matrices, rh 
is an average squark mass, A, B = L,R and i,j are the generation indices. 
The LR(RL) mass insertion can enhance the Wilson coefficients CV 7 and C$ g 
by a factor of nig/nib compared to the standard model contribution. This 
leads to a strong limit of order O(10 -2 ) on the LR(RL) insertions | $lr(rl) I 
from the B(B — > Xsj) j.'UI l3"!2] while the limit on the LL and RR ones is 
rather mild |31II32| . Thus, although larger values for LL and RR mixings are 
allowed, when one considers B — > (j)K, the effect of their mixings are only 
significant in the parameter space where the squark and gluino masses are 
at the edge of their experimental constraints |17|. Motivated by this fact, we 
only concentrate on LR(RL) down type squark mixing in hereafter. Thus, the 
new physics effect is very sensitive to $lr(rl)- 
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In general, these contributions LR(RL) can generate gluonic dipole in- 
teractions with the same as well as opposite chiral structure as the standard 
model. In our analysis we will consider each of them separately. Furthermore 
we will only consider the gluonic dipole moment operator, which is the dom- 
inant operator for this process. 

The effective Wilson coefficient for C% g USY obtained in the mass insertion 
approximation is given by for the same chiral structure as the standard model 



r SUSY ( \ _ u ' rbs <"•; f . , ,.>,,, 

with 



G(x) = — -7 [22 - 20a; - 2x 2 + 16x ln(x) - x 2 ln(x) + 9 ln(x)l , (21) 

3(1 — x) 4 

where x = rn^jm 2 - is the ratio of the gluino and squark mass. 

Using the renormalization group equation one can evolve the coefficient 
Cg g USY from the high scale rriq to the scale rrib relevant for B — ► c/>K decay 

Ci g USY (m b )=r,Ci g USY (m,), (22) 

with 

r] = {a s {mq)/a s {mt)) 2/21 (a s (m t )/a s (m b )) 2/23 (23) 

One can obtain C% g USY for opposite chirality, by adding one more operator 
similar to Os g with (1 + 75) — > (1 — 75) and 8 b ^ R — > S RL . However, in B — > <pK 
process, both LR and RL contribute with the same sign because B and K 
parity are both CP, and the process is parity conserving. 

The effective Wilson coefficient CfJ is defined as CfJ = C 8g + Cjf sy . 

This effective Cg g will contribute to the amplitude M.{B — > (f>K) through 
the function of Equation 1161 C$ g depends on the magnitude and phase 
of the vame °f squark mass (m^) and the ratio x (= m~/m~). The 

variation of C% with x is determined by the function | G{x) | as shown 
in Figure n From this Figure, it is clear that SUSY gluino contribution to 
B — > (f>K first increases with increase in x, and then after some value of 
x = 0.5, it starts decreasing asymptotically with further increase in x. 

The different input parameters and their values used in numerical calcu- 
lation of branching ratio and CP asymmetries are given in Ref. [21] • 

5.1 LR(RL) mixing 

In this section we study the effect of LR(RL) mixing inB-> <j>K process. This 
LR(RL) mixing of the down type squark sector can also affect the B — ► jX s 
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process and B s — B s mixing. Hence we need to take into account the limit 
on LR(RL) mixing parameter $l R ^ R l^ from the above two experimental data 
in the present analysis. In the first case, it has been shown in Ref. [31] that 
from the measurement of B(B — » jX s ) one gets | S^ R , RI s |< 1.0 x 10~ 2 and 
3.0 x 10~ 2 for x = 0.3 and 4 respectively, with rriq — 500 GeV. It is interesting 
to note that the lower the x value stronger the limit on | 5 h £ R ^ RL ^ | , which can 
be explained by the x dependent behavior of the C^ SY . 

The current experimental data on B s — B s mixing is AM S > 14.4 ps _1 ( at 
95% C.L.) We have found that the LR(RL) mixing does not change the 
value of AM S significantly from the standard model prediction in the allowed 
range of | S b L s R , RL) \. 

In our analysis we consider rriq = 500 GeV and take two values of x = 0.3 
and 4.0, which will determine the gluino masses. In Figure|21we show the la 
allowed region in p—ip plane from B — > <pK data on S^k, C^k and B. The gray 
band indicate the parameter space which is allowed by S$k- The area outside 
the two dotted contours is allowed by C^k , while the area enclosed by the 
solid curves is allowed by the B(B — » 4>Ks) measurement. The region (marked 
by Z) in gray band enclosed by the solid curves is the only parameter space 
left in p — 1\) plane which is allowed by the experimentally measured S^k, C^k 
and B within la. 

The Figures |2Ja) and (b) , correspond to contour plots for x — 0.3 and 4.0 
respectively at the scale p = rrib- For x = 0.3, we get two allowed regions each 
at positive and negative values of the new phase tj). On the other hand for 
x = 4.0, we get only one allowed region which lies at the negative value of if) 
and at much higher value of p > 2.2 x 10~ 2 . We have noticed before that the 
constraint on LR(RL) mixing parameter from the B(B — > X s ^/) is stronger 
at x = 0.3 compared to the limit at x — 4.0. This behavior is also reflected 
in the B — > <f>K process, where we find that, for x — 4.0, the la constraint 
from S(j>K, C$k and B is much weaker compared to the constraint shown in 
FigureEIa) correspond to x = 0.3. 

Similar allowed regions are shown in Figure Hfc) for a different choice of 
the QCD scale p — to&/2. One can see that the allowed parameter space does 
depend on p. In this case, both the allowed regions are confined at the positive 
value of ip. For x — 4.0 ( Figure GI<f)), there are no allowed regions. From the 
S^k and branching ratio contour one can sec that the allowed region from 
B — > <f)K require some higher value of p which lies beyond B — > X s j limit. 

Before we conclude this section, we would like to compare our predictions 
with some of the existing literatures onB^ <j>K process [TU1ITT1H7] . 

We agree qualitatively with the results of Ref.^DEj m places where we 
overlap. Similar to our approach, both of these analyses were based upon the 
QCD improved factorization scheme. However, there are some quantitative 
differences between these papers and our analysis. For example, we differ in 
the choice of squark and gluino masses, the authors of the above two papers 
considered degenerate squark and gluino masses, whereas we have considered 
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non-degenerate squark-gluino masses. We have fixed the squark mass at 500 
GeV and considered two values of the gluino masses, determined by the pa- 
rameter x defined earlier. Secondly, we have performed our analysis for two 
values the QCD scale, fi — rrib/2 and vtib- Our results depend strongly on the 
choice of the ratio x and also on the scale p,. However, in a broad sense, we do 
agree that to satisfy B — > <j>K data, one requires | ^lr(RL) l~ 10~ 3 — 10~ 2 . 

In Ref.|10|. authors made a detailed investigation of a scenario in which 
the LR and RR operators co-exist. Moreover, because of the large mixing, 
the calculation was done in the mass eigenbasis with more model dependence 
than ours. It has been shown in this analysis that RR insertion (which arises 
due to a large mixing between sr and 6r) could show sizable effect on S^k, 
but only for very light gluino mass, near the experimental bound. Such a 
large RR mixing also modify AM S significantly which can be observed at the 
Tevatron Run II. In their second case, they have the combination of both large 
right-right and left-right squark mixing (LR + RR) . In this case the squark 
and gluinos could be sufficiently heavy to have no significant enhancement of 
the AM S . 

From our analysis we observe that SUSY leads to a comprehensive un- 
derstanding of B — > (j)K~s data though in a very limited parameter space of 
"LR(RL) ■ ^ n res ^ °f ^ ne P a P er we now explore the consequence of such LR(RL) 
mixing of squarks in the B — ► (f>K* process. 



6B^ 4>K* decay 

In this section we will study the effect of LR(RL) mixing of down type squarks 
to B — > <j)K* process through the gluonic dipole moment operator Cg g . We 
will study the B — > (f>K* process by using the QCD improved factorization. 
Using this method one can compute nonfactorizable corrections to the above 
process in the heavy quark limit. Recently the B — > VV process has been 
computed using QCD improved factorization method j37) . In rest of analysis 
we will follow Ref. gZj. 

The most general Lorentz invariant decay amplitude for the process B — ► 
VV can be expressed as 

M(B(p B ) -> V 1 (e 1 ,p 1 )V 2 {e 27 p2)) oc e^e™ ag^ + bp Bfl PB» + ice^apPiP^-^) 

where the coefficients c correspond to the p-wave amplitude, and a,b to the 
mixture of s and d wave amplitudes. Using these a,b and c coefficients one 
can construct the three helicity amplitudes: 

#oo = [(m% ~ m Vl ~ m 2 y 2 )a + 2 m \p\ m b\ 

H±± = a =F m B p cm c (25) 
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where p cm is the center of mass momentum of the vector meson in the B rest 
frame and my 1 (mv 2 ) is the mass of the vector meson V\{V2). These helicity 
amplitudes H oa and H±± can be related to the spin amplitudes in the trans- 
verse basis (Aq, A\\ , A±) defined in terms of linear polarization of the vector 
mesons: 

A 
A \\ 
A± 



Hqo 
1 

1 

V2 



(#4 
(#4 



H-) 
H-) 



(26) 



Branching ratio 


Data 


Weighted average 


B+ -f cj>K*+ 


BaBar (l2.1±lg ± 1.5) x 10" 6 
CLEO (10.6±2;|±1;1) x 1(T 6 
Belle (9.4 ±1.1 ±0.7) x 10" 6 


(9.9 ± 1.23) x 1(T 6 


B" - <f,K M 


BaBar (ll-lt^a ± 1-1) x 10" b 
CLEO (11.5+if+if) x 10" 6 
Belle (Wt^l) x 10" 6 


(10.6 ± 1.3) x 10" 6 



Table 1. Experimental data of B — > <f>K* decays from BaBar |IU], CLEO and 
Belle and their weighted average. 



The decay rate can be written as 

r(B - y^ 2 ) = [| # 00 p + 1 #++ | 2 + | | 2 ] (27) 

Neglecting the annihilation contributions (which are expected to be small) to 
B — > -ffoo and H±± are given by: 

H 00 = % ° m y g *)/» J ( m | - m | . _ m » ) {mB + m K , ) A™' (m 2 ) 



V2 2m K - 
4m 2 3 p 2 

?™B ± TTlK* 



-Af K *(ml) 



H±± = ^a"(^*)m /J (ffl B ± mjr.) Af K *(m 2 ) 

T^^^-K)) (28) 

where, a n (4>K*) — a% + a J ± ag - (a" ± + a" )/2. The effective parameters 
at appearing in the helicity amplitudes H 00 and H±± given in Ref.[23|^] and 
other input parameters are given in Ref. j24j . 
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In Table 1, we display the experimentally (BaBar, CLEO and Belle) mea- 
sured branching ratios and the weighted averaged values for the B + — > <j)K* + 
and B° — > <pK*°. The theoretical predictions in the SM for two different form 
factor models, the LCSR and BSW models are given in Ref . [3*7] . 



X 


M) 


23 SUSY (in units of HT e ) 


Acp (in %) 


0.3 


(0.4 x lGT*, -0.5) 


23.37tJ-S (21-761J-S) 


-4.7(-4.4) 


(0.4 x 10"*, -0.7) 


21.50i 4 4 ;S (20.17j££) 


-7.0 (-6.5) 


(0.6 x 10"*, 1.5) 


27.461^° (26.821--) 


17.7 (15.7) 


4.0 


(2.4 x 10-^,-0.5) 


24.33l^6 8 (22.651*;™) 


-4.7 (-4.4) 


(2.6 x HT*, -0.45) 


26.98lgi (25.HI476) 


-4.2 (-3.9) 


(2.8 x 10"*, -0.8) 


25.311S;S (23.871^) 


-8.0 (-7.5) 



Table 2. B(B + -> (/>#*+) and _4 C p(B + -> 0A"* + ) at the QCD scale fx = m b for 
LJ? mass insertion for selected points in the allowed p — tj) space. The numbers in 
the parenthesis correspond to the RL mass insertion. The standard model branching 
ratio corresponding to this scale is (6.18jlx!i5) x 10 -6 . The errors are due to ±10% 
theoretical uncertainties in the calculation. 



6.1 LR(RL) mixing contributions to B — > <pK* 

In this section we will study the effect of LR(RL) mixing in B — > <f>K * process. 
This LR(RL) mass insertion can enhance the Wilson coefficient Cs g by a 
factor of mg/mb compared to the standard model contribution in the same 
way as shown in section 5 for B — > <j)K process. Hence, one need to impose 
the constrain on LR(RL) mixing from experimentally measured B(B — > X s j) 
and also from S$k, C^k and B(B — > <pK) as obtained section 5. 

In this scenario, the new weak phase if>, (the phase of the LR(RL) mixing 
) will contribute to direct CP-violating asymmetry Acp defined as : 

= r(B+ ^ <j>K*+) - r(B- ^ <j>K*~) 
CP r{B+ -><f>K*+) + r(B- ^<f>K*-) [ ' 

in terms of partial widths. Recently BaBar and Belle Collaboration has pre- 
sented their measurement of CP violating asymmetries for B° — > (f>K*° and 
B ± -> 0A"*± HOI Hg 

A C p(B° -> (j>K*°) = 0.04 ±0.12 ±0.02, 0.07 ± 0.15±g;g| (30) 
A C p(B ± -> cjjK^) = +0.16 ±0.17 ±0.04, — 0-13 ± 0.29±g;Jf (31) 

where, in each asymmetry result, the first number correspond to the BaBar 
data while the second one correspond to Belle measurement. The standard 
model value for this asymmetry is less than 1%. The new physics (SUSY) 
contributions from the new penguin operator appeared due to LR(RL) mixing 
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B SUSY (in units of 10~ b ) 


Acp (in %) 


0.3 


(0.55 x nr*, i.8) 


(32.45+^) 


19.39+Toi (16.17± ;^ 7 ) 


(0.82 x 10"*, 2.8) 


14.501^ (21.621*;**) 


20.2 (10.961--) 


(0.82 x 10-', 2.9) 


12.39 ;: (19.82^;£) 


15.73 (8.041--) 



Table 3. B(B+ -> 0^*+) and A C p{B + -> 0/C+) at the QCD scale p = m b /2 for 
LR mass insertion for selected points in the allowed p — ijj space. The numbers in 
the parenthesis correspond to the RL mass insertion. The standard model branching 
ratio corresponding to this scale is (14.92i2'°|) x 10 -6 - The errors are due to ±10% 
theoretical uncertainties in the calculation. 



(^ R r RL \) can modify the sign and magnitude of Acp(B ± — ► (pK*^) within 
the allowed parameter space of 5^ R , RL y 

To get the numerical values of B(B+ -> cpK*+), and A C p{B+ -► 4>K*+), 
we fix x = 0.3 and 4.0. Then for a given QCD scale p,, we select some points 
in the allowed parameter space of p — ip plane (as marked by Z in Figure [2J 
for both values of x. In this computation, we include (±10%) theoretical un- 
certainties. 

In Table 2, we present the branching ratio B(B + — > <fiK* + ) and the CP rate 
asymmetry A(B + — > <fiK* + ) for p = m& and selected values of x = 0.3 and 4.0 
for values of p and ip allowed by B — * data for Li? mass insertion (i?L 
is shown in the parenthesis). The branching ratio with SUSY turn out to be 
much higher than the standard model value of 6.181^^5, which is lower than 
the experimental data (Table 1). Even the lower range of theory prediction is 
much higher than the upper range of experimental data within la. The rate 
asymmetry has much less error and is consistently within the range ~ —4% 
to - 18%. 

Similarly in Table 3, we show B and Acp calculated for QCD scale \i = 
mb/2. In this case, there are two allowed regions from the combined B — > <j)K 
and B — > X s j constraints corresponding to x = 0.3. For x = 4.0, there are no 
allowed regions from B — > <f>K data. The standard model branching ratio is 
much larger compared to the one computed at fi = m&. In SUSY, apart from 
the QCD scale p, the branching ratio also depend on the values of p and ip. 
Moreover, the selected points in p — ip plane are different in the two cases. In 
the case, with p = mf,/2, and at p = 0.82 x 10~ 2 and ip = 2.8,2.9 radians, 
with LR mixing, lower ranges of the theory predictions are consistent with 
the upper range of experimental data at one sigma. For the other value of 
p and ip, the theoretical prediction for branching ratio is much higher than 
the standard model theory as well as experimental data. With RL mixing, 
the predicted branching ratio is much larger compared to both the standard 
model prediction and experimental data. The asymmetries for both LR and 
RL mixing case are always positive with less errors. 

We conclude that for some selected points in p — ip plane allowed by B — > 
4>K and B — > X s j at p, = m^/2 provide a satisfactory understanding of 
B — > <pK* process. We also note that, at p, = m& the SUSY contribution 
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to the branching ratio of B — > <f>K* is too large to be consistent with the 
experimental data. 

We have also studied other CP violating asymmetries that can arise in 
vector- vector final state. The set of observables are defined in terms of Ao, A\\ 
and A± as follows |26| . 



Ay 



|A A | a + |A> 



S 



A A 



2 2 
A_u = — Im (A±A* — A_\_A*) , A|| =Re(A\\A% + A\\A* ) 
S u = -Im (Aj_A* + A±A*i) , 27, | = Re (A\\A^ - A||-A*o) 

'<1 



= Rc(J [AlAi + AtAxf) 



PA 



= Im 1A*: A 



V 



p ll0 = -lm( - A^Ao + A* A {1 



pa = -Im -A* A,, 
\P 



(32) 



where A = {0, ||, _L} and the observables where i — {0, ||}, We restrict our- 
selves to the study of helicity dependent CP asymmetry denned as £\\/A\\ 
|2(i|. For the purpose of illustration we select last two sample points from the 
Table 3. At these values of p and ip, with LR mass insertion, the lower range 
of the S SUSY is consistent with the upper range of the experimental data on 
B(B — ► 4>K*) at one sigma. We then compute S\\/A\\ for each values of A 
for these two sets of p and ip and is shown in Table 4. As before, in this case 
also we include ±10% theoretical uncertainties in our calculation. We only 
show the helicity dependent asymmetries for LR mass insertion, since with 
RL mass insertion the SUSY contribution to the branching ratio is too large 
to be consistent with the data. 



X 




■Soo / Aoo 


^ii Mil 




0.3 


(0.82 x 10 -2 , 2.8) 


0.19 ±0.00 


61 +,U)UU 


u -°' -0.013 


(0.82 x HT*, 2.9) 


0.15 ±0.00 


n 71 +(1.1)21) 
u - ' ± -0.026 


U - OJ -0.027 



Table 4. Helicity dependent CP asymmetry at the QCD scale p = rrib/2 for LR 
mass insertion for selected points in the allowed p — ip space. The errors consist of 
±10% theoretical uncertainties. 



7 Conclusions 

In this talk, we considered the SUSY contribution to the gluonic dipole mo- 
ment operator to B — > <pKs process. We found that the LR(RL) mass in- 
sertion can enhance the gluonic dipole moment operator significantly. We 
then used the experimentally measured quantities, such as S^KjC^k and 



New Source of CP violation in B physics ? 13 

B(B — » (f>Ks) to constrain the parameter space of LR(RL) mixing. Interest- 
ingly, we find that the constraints from B — > <f)K data is consistent with the 
B — > X s 7 limit. It turned out that the same enhancement of gluonic dipole 
moment operator could also affect other penguin dominated process, such 
as B — > 4>K* , which is a pure penguin process like B — > 4>Ks- In standard 
model, the predicted Acp(B — > is less than 1%. We calculated such 

asymmetries and also the branching ratio for the set of parameters allowed by 
B — > (f>K data. At /j, — nib, for both LR and RL mass insertion, we observed 
that the predicted branching ratio is well above the experimentally measured 
one. On the other hand, at the QCD scale fi = rrib/2 with LR mass insertion, 
we found that the theoretically computed branching ratio is consistent with 
the data with in one sigma error. At this second choice of allowed parameter 
space of 8 b jS R , RL y we found Acp(B + — ► <pK* + ) in the range 15% to 20%, 
which is significantly higher than the standard model prediction but is still 
consistent with the present data. Finally, we also presented helicity dependent 
CP asymmetries in the same parameter space of S£ R . 
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Fig. 1. Variation of | G(x) j with x{= m|/m?). 
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Fig. 2. Contour plots of S^k, C^k and B(B — > 4>Ks) in p — ip plane for two values 
of x — 0.3 (a, c) and 4.0 (b, d) for LR(RL) mixing with m, = 500 GeV. The scale 
H — m b for Figures (a) and (b), while it is m b /2 for Figures (c) and (d). The la- 
allowed regions of S</,k, B and C^k are two gray bands, area within the solid curves 
and area outside the two dotted contours respectively. 



